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Abstract 



In this paper a super symmetric version of a generalized unitarity cut method 
in application to MHV and NMHV for form factors of operators from the J\f = 4 
SYM stress-tensor current supermultiplet T at one loop is discussed. The explicit 
answers for 3 and 4 point NMHV form factors at tree and one loop level are obtained. 
The general structure of n-point NMHV form factor at one loop is discussed as well 
as the relation between form factor with super momentum equal to zero and the 
_i ■ logarithmic derivative of the super amplitude with respect to the coupling constant. 
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1 Introduction 

Much attention in the past decade has been paid to the study of the scattering amplitudes 
(the ^-matrix) in four dimensional gauge theories, especially in the planar limit of M = 4 
SYM theory. 

It is believed that the hidden symmetries of the M = 4 SYM theory which are respon- 
sible for its integrability properties will completely fix the structure of the amplitudes. 
The hints that the S'-matrix for the M = 4 SYM theory can be fixed by some underlying 
integrable structure were found at weak [U E[ E] and strong [21 0] coupling regimes. 

There is another class of objects of interest in the TV = 4 SYM theory which resemble 
the amplitudes - the form factors which are the matrix elements of the form 

(p*V..,pJ»|0|o>, (i.i) 

where O is some gauge invariant operator which acts on the vacuum and produces some 
state \pi , ■ ■ ■ ,Pn n ) with momenta pi,..-,p n and helicities Ai,...,ArJj. The S'-matrix 



'■Note that scattering amplitudes in "all ingoing" notation can schematically be written as 

(^V..,pMo). 



operator is assumed in both cases. One can think about this object as an amplitude 
of the proses where classical current coupled through gauge invariant operator produces 
quantum state. The example of such process is 7* — > Jet's in perturbative QCD [5] (see 
also [6] |7] for recent results) where we take into account all orders in a s but the first order 
in a em . The amplitude of such process is given by the matrix element of the following 
form: 

(^•••pMj^Io), 

where j2% D is the QCD quark electromagnetic current operator. 

The two-point form factors in M = 4 SYM were studied long time ago in [8] and 
recently in [10]. Using the M = 3 superfield formalism the form factors of none gauge 
invariant operators (off-shell currents) at tree level were derived in [9]. Recently, the strong 
coupling limit of form factors has been studied in [IT] and the weak coupling regime in 
[T2"| [T3"l Q3J [T51 [T5]. Also different regularizations for form factors were discussed in [17]. 
The motivations for the systematic study of form factors in M = 4 SYM are 

• it might help in understanding of the symmetry properties of the amplitudes [H [2]. 
It is believed that the symmetries completely fix the amplitudes of the J\f = 4 SYM 
theory and it is interesting to see whether they fix/restrict the form factors as well; 

• the form factors are the intermediate objects between the fully on-shell quantities 
such as the amplitudes and the fully off-shell quantities such as the correlation func- 
tions (which are one of the central objects in AdS/CFT). Since the powerful com- 
putational methods have appeared recently for the amplitudes in M = 4 SYM ( see, 
for example, [HI [19] and [20j|2T]), it would be desirable to have some analog of them 
for the correlation functions [22]. The understanding of the structure of form factors 
and the development of computational methods might shed light on the correlation 
functions; 

• also, it might be useful for understanding of the relation between the conventional 
description of the gauge theory in terms of local operators and its (possible) descrip- 
tion in terms of Wilson loops. The latter fact is the so-called amplitude/Wilson loop 
duality which originated for the case of M = 4 SYM in [231 EH [25] . This duality 
was intensively studied in the weak and strong coupling regimes and tested in dif- 
ferent cases, and its generalizations to the non-MHV amplitudes were proposed in 
[25] . Moreover such dual description for amplitudes of M = 4 SYM together with 
the developments of OPE technique for Wilson loops [27J [28] led to the proposal of 
the equation [29] which in principle should define the whole M = 4 SYM S-matrix 
for any value of coupling constant. Note that similar equation for the amplitudes 
may be derived from "twistor space" (see for example [30]) point of view [31]. It is 
interesting to investigate whether such dual description for form factors in M = 4 
SYM exists and if it is possible to formulate this equation for form factors. 

To make progress in the above-mentioned directions, the perturbative computations 
at several first orders of perturbative theory are likely required. 



The aim of this paper is to continue investigations of the form factors jH] of the 
operators from stress tensor supermultiplet in M = 4 SYM. We will use the developed in 
[13] formulation of form factors in M = 4 on-shell momentum superspace, which allows 
us to consider form factors with different types of particles in (p 1 1 , . . . , p^ n | external state 
in M = 4 covariant manner. We will use the generalized unitarity technique to study the 
structure of NMHV sector at one loop. First we are going to discuss how the generalized 
unitarity technique works for form factors for MHV sector at one loop in M = 4 on-shell 
momentum superspace. Then we will continue with the NMHV sector. We will perform 
explicit computations of 3 and 4 point NMHV form factors at one loop and will discuss 
the structure of the general n point situation. We make a brief comment on the relation 
between form factors with operator insertion with zero momentum and amplitudes. 

2 Amplitudes and form factors in on-shell momen- 
tum superspace 

2.1 Super form factors of chiral truncation of M — 4 stress tensor 
supermultiplet 

To describe the M = 4 SYM stress tensor supermultiplet it is convenient to use standard 
N = 4 coordinate superspace 

J\f = 4 coordinate superspace = {x aa , 9 A , OAa}, (2.2) 

where x a6i are bosonic coordinates and 0's, which are SU(A)r vectors and Lorentz SL(2, C) 
spinors, are fermionic ones. The Af = 4 supermultiplet of fields (containing (p AB scalars 
(anti-symmetric in 5 , L r (4)^ indices AB), i[) A , ^) A fermions and F^ v - the gauge field strength 
tensor, all in the adjoint representation of SU(N C ) gauge group) is realized in M = 4 co- 
ordinate superspace as a constrained superfield W AB (x,9,9) with the lowest component 
W AB (x, 0, 0) = cf) AB (x). W AB in general is not a chiral object and satisfies several con- 
straints [321 E] : a self-duality constraint 

W AB (x,6,6)=W^(x,6,6) = ±e ABCD W CD (x,6,6), (2.3) 

which implies (f) AB = 4>ab — \t ABCD 4>CD and two additional constraint^ 

D%W AB (x,e,e) = -p lA D a L W B ^ L (x,d,6), 

D^ c W A)B {x,6,6) = 0, (2.4) 



,*] denotes antisymmetrization in indices, while (*,*) denotes symmetrization in indices. 



where D A is a standard coordinate superspace derivative^. Note that in this formulation 
the full M = 4 supermultiplet of fields is on-shell in the sense that the algebra of the gen- 
erators Q A , Qbo of the supersymmetric transformation of the fields in this supermultiplet 
is closed only if the fields obey their equations of motion. The J\f = 4 SYM stress tensor 
supermultiplet T AB is given then by 



^AB 



Tr (W AB W AB ) . (2.5) 



We will consider in this article the chiral truncation of the M = 4 SYM stress tensor 
supermultiplet (which contains only self-dual part of full multiplet) rather then the super- 
multiplet itself [H] . The main reason for this is that the chiral truncation has the off-shell 
description in terms of superfields on M = 4 superspace i.e. the component fields in such 
truncated multiplet are arbitrary and the chiral part of the algebra of supersymmetric 
transformations of the component fields can be still closed without any constraints on the 
component fields. Note that the off-shell description for the full Af = 4 supermultiplet in 
any superspace is unknown. 

To describe this truncated supermultiplet one has to break SU(4:) R group into two 
SU(2) and £7(1) 

SU{A) R -)• SU{2) x SU{2)' x 17(1), (2.6) 

so that the index A of _R-symmetry group SU{A)r splits into 

A -» (+o|-o'), (2.7) 

where +a and —a' correspond to two copies of SU(2) and ± correspond to the U(l) 
charge. We will not write the U(l) factor explicitly hereafter, and will use the notation 

(+a\-a') = (a\d). (2.8) 

After that one has to take the particular (ah) projection of W AB that depends on half 
of the Grassmann coordinates [H] (this can be seen from the (12.41) ): W (x, 6 C , 6c)- The 
truncated stress tensor supermultiplet is then given by 

T ab = Tr (W ab W ab ) . (2.9) 

v ' 0=0 

To describe external states in M = 4 covariant manner it is convenient to use M = 4 on- 
shell momentum superspace [53]. This superspace is parameterized in terms of SL(2, C) 
spinors A a ,Aa,a,d = 1,2 and Grassmannian coordinates r] A ,A = 1,...,4 which are 
Lorentz scalar s and SU(4)r vectors 

On-shell M = 4 momentum superspace = {A Q , Aq, r/ A }. (2-10) 



3 which is D A = d/d9% + i6 A6l d/dx a 



In this superspace the creation/annihilation operators 

{g-, T A , <P AB , T A , g+}, 

of the M = 4 supermultiplet, for the on-shell states which are two physical polarizations 
of gluons \g~), \g + ), four fermions \T A ) with positive and four fermions \T A ) with negative 
helicity, and three complex scalars \4> AB ) (anti-symmetric in SU(4)r indices AB ) can be 
combined together into one M = 4 invariant superstate ("superwave-function") \Qi) = 
^|0>: 

la) = U + V A r hA + ^V A V B <Pi,AB + y% B V C eABCnT? + ^V A V B V C V D e A BCD9-) |0>, 



(2.11^ 



where i corresponds to the on-shell momentum p l a6l = A^A^, pf = carried by the particle. 
The n particle external state \Q n ) is then given by \Q n ) = YYt=i^i\fy- 

The form factor JF n of the truncated stress tensor supermultiplet for general n particle 
external state is then given by: 

F n ({\,~\, V },q,d a ) = (Q n \T ab (x,6 a )\0), (2.12) 

where {A, A, 77} is short notation for (Ai, Ai, 771 . . . A n , X n , r) n ). Here we are considering 
colour ordered object JF n . The physical form factor T^ 8 ' in the planar limito should be 
obtained from J 1 ^ as: 

H h "'({A, A, v], q, 6 a ) = (2vr) V" 2 2 n/2 £ Tr(t a ^ . . . t a ^))JF n {o{{\, A, V }), q, 9 a ), 

(2.13) 

where the sum runs over all possible none-cyclic permutations a of the set {A, X,f]} and 
the trace involves SU(N C ) generators t a in the fundamental representations. The normal- 
ization Tr(t a t b ) = 1/2 is used. 

Performing Fourier transform for bosonic coordinate x aa — > q a a and taking into ac- 
count that J- n is chiral, translationally invariant and T ab is 1/2-BPS we see that JF n should 
satisfy the following conditions [14] : 

where generators of supersymmetry algebra (P aa , Q a a -, Q a a , Qaa, Qaa) acting on T n are given 
by 



4 translations P aa = — 2_, KxKt + 1o 



i=X 



l g — > and N c —y 00 of SU(N C ) gauge group so that A = g 2 N c =fixed. 



n d 

4 supercharges Q a a = ~ E A «^ + ^ 



dO C a 
1=1 " 

n d 

4 supercharges Q a a = - E,W + — , 

t=l o 

4 conjugated supercharges Q ad = - ) A^-— + 0"g o 



n ~. (9 
4 conjugated supercharges Q ad = -J^ ^^i + ^?9ad- ( 2 -!5) 

i=i 
This relations imply that T n takes the following form [14] : 



! 9r/| 



J-„({A,A^},?,e B ) = ^(^A^-g^e^^^^^A,^}^ 



8=1 



^ = A , W + A' r ( 4 ) + ... + A'i 4n - 8 ), (2.16) 

where 

n n 2 

«S = E « > 9S = E X ^t ^h(«S /4 ) = II ^V^' ( 2 - 17 ) 

i=1 i=1 0/0,6/6=1 

and Af„ m ' are the homogenous SU(4)r invariant (more accurately SU(2) x SU{2)' x t/(l) 
invariant, they can be written in SU(4)r covariant form in M = 4 harmonic superspace 
formulation) polynomials of the order of 4m. 

Assigning helicity A = +1 to |fl$) and A = +1/2 to 7/ and A = —1/2 to 0^, one sees 
that J- n has an overall helicity As = n, 5q R has As = 2, exponential factor has As = 
so that X n has As = n — 2, X n has As = n — 4, etc. A*„ , A™ etc. are understood 
as analogs [33] of the MHV, NMHV etc. parts of superamplitude i.e. part of super form 
factor proportional to Af„ will contain component form factors with overall helicity n — 2 
which we will call MHV form factors, part of super form factor proportional to Xn will 
contain component form factors with overall helicity n — 4 which we will call NMHV etc. 
up to X n n ~ overall helicity 2 — n which we will call MHV. 

It is convenient to perform transformation from Q a a to the set of axillary variables 

{x' a ,v a ,Cv" a y 

r 2 

f [...]= /^exp(^E«)[-4 (2.18) 

After such transformation we can write T[J-" n ] as (let's use the notation X a r] a + X a r] a = 

n 

Z n ({\,\,v},q,{la}) = f[7 n ] = 5\Ys X ^-q a M R (<l a a + lMR{(l a «) Xn- (2.19) 

i=\ 

7 



The algorithm of obtaining component form factors from this supersymmetric expression 
was discussed in [TJ]. 

Let's make a comment about total U(l) charge of T n and TfJ-^]. The Tab operator 
carries (—4) charge, the (Q n \ carries charge (0) so the J- n form factor should carry (—4) 
charge, which indeed true and can be seen from fl 2 . X 6 j) : Sq R {cfy carries (—4) charge, 

while e e °- q °' and Xn are neutral. The T transformation (the integration measure d 4 0^) 
also carries (+4) charge, so that Z n = T[J- n ] is neutral with respect to U(l). Note that 
this will be no longer true for the form factors of operators from different supermultiplets. 
The formulation of form factors discussed so far lacks of explicit SU(4)r covariance. 
SU(4)r covariance can be restored in the M = 4 harmonic superspace formulation. How- 
ever such formulation does not give us any computational benefits for the purpose of 
our computation and all results obtained in our none covariant formulation can be easily 
translated to SU(4)r covariant formulation. We will discuss such formulation briefly in 
appendix. 

2.2 Generalized unitarity for form factors at one loop 

In general at one loodj level Zn can be decomposed as combination of all possible scalar 
boxes (B im ,B 3m ,B 2mh ,B 2me ,B lm ) and triangles (T 3m ,T 2m ,T lm ) integrals : 

y(l) _ \ y /~t&rn r>Am _i_ /i3m r?3m _i_ /i2m/i r>2mh < /~i2me r>2me _i_ /ilm nlm 

lb / _J % % % % L % it It L % 

i 

+ J2 C 3 mT 3 m + C j mT j m + C j mT j m+ P erm -> ( 2 ' 2 °) 

3 

where the sum runs over all possible distributions of the ordered set (pi, . . . , p n ) of individ- 
ual momenta between vertexes of the scalar integrals, while the position of the momentum 
q carried by the operator is fixed, perm, corresponds to the cyclic permutations of the 
(pi, . . . ,p n ) set of the momenta of external particles. The latter is necessary due to the 
fact that while we are considering object that is colour ordered in the colour space of 
external particles the operator is colour singlet hence the momentum q carried by the 
operator can be incepted at any position in the colour ordering [13] . This is equivalent to 
the consideration of all possible permutations of external momenta, while the position of 
the operator momenta q is fixed. In general we can write scalar box integral as: 

f d D l 1 

B K lKlKlKt - J {2n)D [2{Ki - l)2{Ki - ^ - /)2(/ _ Ki)2 , (2.21) 

where J2i=i K&a = Yh=i KxKt~1aa- The particular scalar box integrals B 4m , B 3m , B 2mh , 
B 2me , B lm are defined then as: for Bt% „ 2 „ 2 „ 2 all K 2 + 0, for B*% k2 w2 K 2 = 0, for 



^Hereafter we do not write common one loop factor (\in D / 2 rr)/(2Tr) D explicitly. See appendix. 



B*™* K{ = K l 2 = 0, for Bpi e K2 K( = K* = 0, for BVS K( = K\ = K l z = 0. For 

3' 4 2 1 4 4 

triangle integrals we use similar notations: 



TkIkIk! - J {2 n)Dp(K 1 + iy(l-K 3 r (2 ' 22) 

where £? =1 K^ = ££=i KK - fob for T*?^ # all Kf ± 0, for T*?^ if? = 0, for I*J 

if? = #2 = 

The dependence on the helicities of the external particles as well as the type of operator 
are encoded in the C^ coefficients. The C& coefficients are Grassmann polynomials and 
in general, as was explaned earlier, should have the form: 

Cu ~ S* GR (£ + r a ) 5 4 cn (C) {CT + C^ + ...+ Cf~ 8) ) , (2.23) 

where C« are the homogenous SU(4)r invariant polynomials of the order of Am. For 
example coefficients before scalar integrals for NMHV form factors will be proportional to 
^gr (o.a + 7a) $gr (la) Q ■ ^he analytical answers for all types of one loop triangles and 
boxes are known (see [31] for review) and therefore the problem of computation of Zn 
reduces to the determination of Ck coefficients. The latter are computed in the unitarity 
based methods by comparing the analytical properties of both sides of the relation (12.201) 
viewed as the functions of Mandelstam kinematical invariants of momenta of external 
particles. 

To obtain the values of the coefficients C, before box integrals it is very convenient to 
consider quadruple cuts [SSI EH]- Such cuts are unique to each box integrals and stop the 
loop momenta flow. Therefore the quadruple cut completely determines the value of the 
Ci coefficient for the chosen box integral. One can schematically write that 

1 r 4 

Q = -J2 Jl d2 vtd 2 vt Z tree x A tree x A tree x A tree , (2.24) 

±5 •* i 

where Yl±s corresponds to the summation over the solutions of the on-shell and momen- 
tum conservation conditions [361 EI]- We will not write ^2± s sums explicitly in the most 
cases in the next chapters. Note that as in the case of amplitudes [36] for the MHV and 
NMHV form factors we will not need the explicit form of the solution. Also note that 
the summation over the states (types of particles) that runs through the cuts is "hidden" 
in the Grassmann integration ("supersums"). Since we are working with M = 4 SYM at 
one loop all cuts can be evaluated in D = 4 with O(e) accuracy. Z n and A n correspond to 
form factors and amplitudes stripped from the overall delta function 5 4 (J^™ =1 A^A^ — q a a) 
of momentum conservation. The exact types of form factors and amplitudes entering the 



6 Each scalar box can be uniquely specified by its leading singularity. The latter are obtained by 
cutting all four scalar propagators in the integral |35) . 



expression are determined by what type of one-loop form factor (MHV, NMHV, etc.) and 
what particular cut we are considering. 

The situation with the coefficients before triangle scalar integrals is a little more in- 
volved. Though triple cuts are unique to each triangle integral the triple cut do not stop 
the flow of the loop momenta completely and leaves one parameter integral j dt. However 
one can construct algorithm that allows one to fix the corresponding coefficient using the 
form of the triple cut intergrand [3"T| HS] . One can parametrise the momenta of particles 
l? a which crosses the cuts and associated spinors A^, A^ in terms of combinations of exter- 
nal momenta and the t parameter which is the remainder of loop integral (see appendix 
and [37] for details). Then one can obtain the following relation 

Cj = Inf t [^ J2 [ II A? A? (Z tree x A tree x A tree )}{t) _ , (2.25) 

±s J i t=0 

where Ini t means that one takes expansion in t at t — > oo and separate term proportional 
to t°. In other words the coefficient before scalar triangle integral is given by the first term 
in the series expansion of the corresponding triple cut integrand in t at infinity [37J EH] . 
In the case of MHV form factors we will not need the explicit form of the solutions ±S, 
while the case of NMHV form factor is more involved. We will use IR properties of the 
form factors (see below) in some cases to adjust the value of the coefficients before triangle 
integrals instead of the direct computations. 

2.3 Grassmann delta functions 

Throughout this article we will need different types of the Grassmann valued delta func- 
tions. The following notations will be usecO for 8q R type delta functions introduced 
earlier: 

^(Xi) = Sh R (X^ R (Xt), (2.26) 

while for general M we have: 



M 



5 2M {X*) = \\ e^X^Xl (2.27) 



A,B=1 

We also will use the Grassmann delta functions of another type: 



^(£^a) = II(XyC0, (2.28) 

i A=l i 



7 The relevant for us cases are A/" = 2, A/" = 4. 
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where A index runs from 1 to Af, r\f are Grassmann variables and Ci are bosonic ones. 
For such delta functions we will use the following notation: 

PCEvtG) = ^(^VtC^iJ^VtQ). (2.29) 

i i i 

We will also use for saving space the notation d 2 r] a d 2 r] a = d Ar q A for the integration measure. 
In computation of the corresponding coefficients before scalar boxes and triangles we 
will be performing multiple Grassmann integrals with Grassmann delta functions. The 
following relation is extremely usefull: for some q A (A runs from 1 to Af) q A = Y^=i KxW 
one can write the following expansion over some basis 

Qa=^^ + ^7rl, l<*<n, l<m<n, m^l, 
{ml) {lm) 

where A;, A m should be linear independent. The latter relation implies that 
*%£) = (lm)"*" U + £ ^rA Z» U + £ H^) ,i±hi±m. (2-30) 



For example using this relation one can immediately show that [TJ 

' d«rj h d% 2 6^ (X l X + X^rj A + Qi) 5 2Af (AL 1 ^ + \ h a < ~ P*) 

= {kl 2 f '5^ (P A + Qi) , (2.31) 

which is important relation for the two particle supersums. 

We also want to note that the the computation of integrals over J d 2 vjl and J d 2 r)°' 
in the quadruple and triple cuts may be different in details, but can be performed in 
such a way that leads to the same bosonic coefficient and slightly different Grassmann 
delta functions. One can formulate the following rule to simplify computations: one takes 
integrals over J cPrjfcPrji formally replacing 

te + rM R (K) -► P&2 + it) (2-32) 

and integrating over f d 4 r]f. After the integration we have to put 7^ = 0. 



2.4 Tree level MHV and MHV amplitudes and form factors 

We will need as the building blocks in the computation of the coefficients before scalar 
integrals in the MHV and NMHV case at one loop level several explicit expressions for 
the tree level form factors [H] and amplitudes [36] : 

r7tree,MHV _ <:A(ST^ \a\a , aa\ "GR Via + la) "GR \1a) /n QQ \ 

Zn -*{^\\+<l ) (12) ; ; ; (nl) • (2-33) 

11 



This expression is valid for any n > 2 without any kinematical constraints on Aj, Aj 
variables. For the MHV amplitudes we have: 

n ~ X 8 (n A \ 

a tree, MHV _ r4/\ ^ \a\a\ u Via.) to QA\ 

An ~ 6 ( 2-, A * Ai) (12>...(nl>- (2 ' 34) 

This expression is valid for any n > 4 without any kinematical constraints. While for 

n = 3 this expression exists only for the complex values of the momentao in (H ) 

signature. This implies the following kinematical constraints on Aj, Aj variables [36] : 

MHV3 : \{ik) = —Xj(jk), [ij] = 0, any i,j,k at the same vertex. (2.35) 

We will also need the expression for the three point MHV3 amplitude [36] : 

.tree, mhv _ ,4/V" \ a \a jHVi[ 2 ^\ + cycl.perm.) 
A - S (L. X i\) [12] [23] [31] • (2 ' 36) 

This expression also exists only for the complex values of the momenta, that implies the 
following constraints: 



MHV3 : Xi[ik] = —Xj[jk] } (ij) = 0, any i, j, k at the same vertex. 

Note that this expression is 4'th degree in r^'s while all MHV amplitudes and form factors 
are of the 8'th degree. We also will need the form of Z 3 ree ' which is NMHV form 

factor at the same time in full analogy with n = 5 point amplitude. We will discuss the 
structure of it in separate section later on. 

The discussed above kinematical constrains on Aj, Aj spinors immediately lead to the 
fact that some configurations of MHV and MHV vertexes give vanishing result [36] . See 
figffl 

3 MHV warm-up 

As a warm-up before computations of NMHV form factors we will discuss how generalized 
unitarity works in MHV case. MHV form factors should be the lowest components in r^'s 
expansion of Z n of Grassmann degree 8. This implies the following configurations of MHV 
and MHV vertexes 

MHV x MHV x MHV 3 x MHV 3 

for the quadruple cut integrand, and 



MHV x MHV x MHV 3 

for the triple cut integrand. Taking into account that configurations of vertexes depicted 
on figfj] vanish we conclude that the only contributing cuts are those depicted on figj2] 
and figEJ They correspond to the B 2me (and B lm as the limiting case) and T 2m scalar 



8 One can make analytical continuation to the real values of the momenta of external particles in the 
final expressions. 
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Figure 1: Vanishing configurations of MHV3 (grey) and MHV3 (white) vertexes. Dark 
grey blob corresponds to other parts of " diagram" . 



integrals. So the MHV form factor should have the form 

z (i),mhv = J2 C 2 ™B 2me + Y, Cf 1 T I m + Verm., 



(3.37) 



(B lm scalar integrals can also appear as the limiting case of B 2me ). Let's start with the 
quadruple cuts for B 2me integral. I.e. we are going to compute the corresponding C 2me 
coefficient before such integral. This computation is very similar to those in [36J. We will 
treat the configurations where the vertex Z 2 ree ' is present and configurations with 

gtree,MHV ^ n > 3 separately The usefulness of such treatment will be clear later on. 
We have for the configuration A) which corresponds to B 2 ™ e g scalar integral (the 

notation 

s-l 



E 



Pi 



Vr. 



.s-l 



(3.38) 



is used): 



S2...s-l,Sl...s 



X 



1 1 A/ 



A 8 4 (Vi[yi} -%j[M] + Vhi^h}) 8 4 (Vs[hh] -Vhfas] +T]i 3 [sh}) 



t=i 



[U2W1W1I] [sk][kh][hs] 



(hh) . . . (nh) 



■vs-l 



X 



(hh) ■ ■ ■ (2h) 



i) 



(3.39) 



9 We will suppress a SL(2,C) indexes indices in arguments of Grassmann delta functions in some 
cases. 
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V, 





B) 



Figure 2: All possible cuts for the box scalar integrals to MHV form factor. Dark grey 
vertex is MHV form factor, grey vertex is MHV amplitude, white vertex is MHV3 ampli- 
tude. 





A) 



Figure 3: All possible cuts for the triangle scalar integrals to MHV form factor. 

Performing Grassmann integration with first two 5 4 and then integrating the remaining 
delta functions with the use of (12.301) and kinematical constraints of the MHV vertexes 
( f)2.32p also can be used to simplify algebraic manipulations) we obtain 



C 



2me 
S2...s-l,Sl...s 



7tree,MHV 



(12)(s-ls)(ss + l)(nl)[l\hh\s] 



[l\hU\s] [l\kh\s] [Ui][Uh] (s-ll 3 )(2l 2 ) [hU] (rdiXs + lh)' 

(3.40) 



This expression can be transformed with the use of momentum conservation conditions 

lematical 

rtree,MHV ■ 



in each vertex and MHV kinematical constraints to the form 

1. 



C 



2me 
S2...s-l,Sl. 



7tree,MHV 



-Trikhhh) = 
\{{h ~ kf{h - h? ~ (h ~ h?{k ~ k) 2 )- (3.41) 
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This result is identical to the MHV amplitude case [36J. Note also that as was in the 
amplitude case the Eis can ^ e evaluated without the use of explicit solutions for l t . So 
we have that the coefficient before B 2me „ scalar integral takes the form: 

i 

/~i2me iytree,MHV _ a 2me 

Ly s 2 ... s -i,s 1 ... s ~ ^n 2 *2...«-l,Sl...« 



2me 

-i,«l.... = s 2...«-l«l...« - »1 

)2rrae 



A X!U<1.... = S 2 ... s _ lSl ... s - 8i.„ a _i82... a . (3.42) 



The case of configuration B) which corresponds to the B 2me 2 scalar integral is slightly 
different in computational details: 

i 



°Sw =j\\ d 






[li 2 ][/ 2 /i][/il] MIWsHfen 



___^ JLl^ + A^-A^) 

(W 2 (hh)...(2l 2 ) 

(3.43) 

but leads to essentially the same final result: 



r 2me 7 tr,.MUY (12) (n - In) (nl)[l\hk\n}^ 



2me _ r?tree,MHV 

* 2 "- 1 ' 92 "" [lKi^ln] [l/ 2 ][/ 2 /i] (n-ll 3 )(l 3 l 2 )(2l 2 ) [hk][hn] (hh) 



(3.44) 

)2me 

, uau UC OlllipilllCU. BU lllctl tllC UUClllUlCllU UC1U1C -L. 

scalar integral takes the form: 



which, as in the previous case, can be simplified so that the coefficient before B 2 



ITV 



/i2me ytree,MH V _ A 2me 

S2... n -i,q 2 ~ n 2 s2 --™- 1 '? 2 

A %".n-r,*> = S 2...n-ig 2 - Sl...,-lS2...n- (3.45) 

Now let's consider triple cuts. Let's consider the A) configuration which corresponds 
to the coefficient before T 2m 2 scalar integral. We have for the corresponding triple cut 
integrand 

M vihk] - y h [lil] + mAlh]) SgrCZi^ ^vt + 7 a fe(Ez 2 ,z 3 W) 
[lfeHWilPxl] (Ws) 2 

(W • • • W ' ( } 

Performing Grassmann integration in the same fashion as in the quadruple cut cases (first 
integrating 5 2 , then with the use of (12.301) and kinematical constraints of the MHV vertex 
integrating 5 8 ) we obtain for the triple cut integrand: 

7 tree,MHV (12) (nl) [1 ^ |l 3 ) 4 

(^ 2 ) 2 [l/ 2 ][/ 2 /i][/il] (nh)(hh)(h2)- { ' } 
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This expression can be further simplified with the use of the momentum conservation 
conditions associated with each vertex to the forni 10 !: 



Hree,MHV Tr ( 2l 2hh) 



(3.48) 



(J»2) ' 

where the trace can be evaluated, with the use of momentum conservation conditions and 
kinematical constraints of the MHV vertex to the form 



Z n 



tree,MHV 



i(2(.,„.- g ') + g ^--^ 1 + ^ ). (3.49) 



Now applying the following parametrization for lf a : lf a = A^ a + tA^ a + l/tA 3 a , where 
A°2 3 are some constants that depend on external momenta (see [37] and appendix for 
details), we see that 

n 2m, _ ytree,MHV^-j„f \c,( „ „ 2 \ _l V^ (^^...ra — (2, 1 + q)q , 

U s 2 ...n,q 2 - Z n -im t [Z{S 2 ...n ~ 9 ) + } J 



4""—" * ' ^ (l 3 2) 



t=o 



= Z^ MHvl -{s 2 ... n -q\ (3.50) 

so that the coefficient before T s 2m 2 triangle scalar integral takes the form: 

i 

^2m ytree,MHV _ A 2m 

s 2 ...„,q 2 n 2 s 2...n,q 2 

Aj.n,^ = «2...n-? 2 - (3.51) 

Now let's consider the 5) configuration which corresponds to the coefficient before scalar 
integral Tf m . We have for the corresponding triple cut integrand 

yr ,4 J\vihh] ~ WU] + %[U 2 ]) ^ 8 (EIUi W + A^ ~ ^<) 
fj ^ [l/ 2 ][yi][/il] (hl 3 )...(nh) 

<%fl(£U ^ry? + 7 Q + A^X ~ A t ,q£)<%fl(EU W + V?| - Vj|) 

Performing the Grassmann integration with delta functions we obtain: 

(3.53) 



HreeMHV (12) (nl) (ss + l)[l\h\k) 4 



(l 2 2)(sl 3 )(l 3 l 2 ) [ll 2 ][l 2 h][hl] (l 3 s + l)(nh)(hl 3 ) 



10 We use the following notations for the scalar products of momenta of external particles (pi...i,Pj...k) 
(i + ...+l,j + ... + k). 
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This expression can be simplified with the use of the momentum conservation conditions 
associated with each vertex to the form 



Z, 



treeMHV l V^ ( Tr{l x l 2 s + l/ 3 ) Tr{hl 2 sl 



n 



1^ Trihhs + lh) Tr{hl 2 sk) \ 



4 ^ V (hs + 1) (hs) 



The traces can be evaluated as in the previous case and we get 



4^U^ + 1) {hs) 



fe=s+l fc=s+l fc=s+l 



(3.55) 



Note that .Dy depends only on the external momenta. So using the parametrization for 
If a : lf a = A" a + tASf* + l/tA% a , where A\\ 3 are some constants that depend on external 
momenta, we see that 



Cl™ S1 = Z^ MHV - A yird t [ -^- + £±] =0. (3.56) 

ss+i...n,si... n n ^/_^ s L(/ aS + l) (; 3S ) J t _ n 



4=0 



±5 

So we conclude that potentially contributing integrals T 2m 8 does not actually appear. 
This fact was first established in [T3] . 

The latter fact may be puzzling. The resolution of this puzzle comes from the observa- 
tion that Z 2 ree ' and Z^ ee ' , n > 2 vertexes that enter A) and B) triple cuts have 
different number of A/ t spinors in denominators, so that the Z 2 Tee ' vertex is in some 
sense singled out (that's why we treated A) and B) cases separately for the quadruple 
cuts) 

ytree,MHV \ r?tree,MHV \ 

~ (W n ~ (ik)(hh)(jhY 

Using explicit solutions for A^ in terms of t and external momenta [37] for T 2m triangles 

one sees that 

1 

t 



1 



r^j 



(hh) 2 

so if we take into account the whole expression we have none vanishing contribution in 
t — > oo limit, while 



1 ~r>. 



(ih){hh)tjh) 

so in the t — > oo limit we get zero. We verified that this is the general pattern for MHV and 
NMHV cases for the T 2m and T 3m integrals, so on general grounds we can conclude that 
the only allowed triangle integrals should necessary contain one massive q 2 leg in MHV 
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and NMHV sectors. This fact reduces the number of necessary triple cuts significantly. 
Note also that the properties of Z 2 ree ' vertex in the t — > oo limit resemble those of 

the z — > oo limit in the BCFW recursion. One also can use only z^ ee ' MHV , n > 2 but 
not Z 2 ree ' in the BCFW recursion [TBI . 

We see now that the contributing type of scalar integrals in the MHV case are 
B 2 ™ e , Sl , T 2m - (and B 2 ™ e „ 2 as the limiting case of 5 2me , ). The coefficients 

S2...S— l,Sl...s' S2...ni1 V S2...n— 1,9 ° *2...s — 1 i*l...s ' 

before these integrals are given correspondingly by f l3.42[) and (13. 50 p . 

Note also that the combinations of scalar products of external momenta A's in the 
coefficients before corresponding scalar integrals (boxes and triangles) match the scalar 
integral in such way that the coefficients before the 1/e 2 IR pole (the use of the dimensional 
regularization is implied) will have the form ~ A _1 /e 2 (see appendix). This allows us to 
define dimensionless functions E % and T\ just as in the case of amplitudes [36], for all 
types of boxes and triangles as the result of evaluation of the corresponding scalar integral 
through 0(e) multiplied by the corresponding A coefficient. 

As an example, for three point MHV form factor Z\ one can obtain, in precise 



agreement with [TBI [TI] : 



y{l),MHV ir?tree,MHV 



B lm (l,2,3\q 2 



+ 
+ 



T 2m (l|g 2 ,(2 + 3) 
T 2m (3|g 2 , (1 + 2) 



♦5» 

2\ 



lm i 



l,3,2|g 2 ) + ii3 
T 2m (2|g 2 , (1 + 3) 2 



lm 



(2,l,3|g 2 



(3.57) 



Here we write the ordering of massive /massless legs explicitly, using the convention: 

(/'(massless legs| (massive legs) 2 ), 

where Q % is dimensionless function based on the type of scalar integral under consideration. 
The IR divergent part of this result is given by: 



/zl 



IR 



1 3 
72 2^ 



Sii+1 



(3.58) 



while finite part is: 



p(l),MHV ir,tree,MHV 



fin 



-Li 2 1 



-Log 
2 6 



«12 



— I Log 

s 2 J 2 



Note that the IR divergences in the sum of B 
bine in such a way that 



2me 
S2...3- 



Li 2 1 - 



■523 



S23 



LU 



— - -Log — 
S31/ 2 V s i2 



S31 



S31 

7T 2 

2 



■l,Sl. 



Z(1),MHV 



IR 



rtree.MHV 

J n 



1 n 
72 2-^ 



(3.59) 
T 2m 2 will always [13J com- 



s ii+l 
T 2 " 



(3.60) 
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Figure 4: Example of the dual contur for the MHV n=3 form factor. 



This is in fact the consequence of the fact that IR poles should cancel in IR finite ob- 
servables [38] such as inclusive cross sections [391 EH El H21 03], energy flow functions 
[4*4"1 l4"5l 14*6] . etc. based on form factors [H H7] and we expect that similar behavior will 
take place for all types (MHV, NMHV, etc) of form factors 



Z™ 



7tree 



IR 



1 " 
72 2.^ 



F 



(3.61) 



It was noticed in [13] that it is likely possible to find the perturbative description of 
the form factor in terms of the periodic Wilson loops. This fact suggests that the use of 
the dual variables [33] as in the amplitude case will be useful though the dual conformal 
properties [T2l [i~4"j of form factors remain obscure. We will use dual variables as compact 
notations in some cases. One introduces dual coordinates x? a as 



- ' .,- 



£<' 



(3.62) 



Note that in the case of the periodic contour within one period there are n labels for the 
momenta of external particles but n + 1 label for the dual Xi points. For example we 
can write the momentum carried by the operator q as q = Y^i=\Vi — x in+i- See figJH 
for the n = 3 case. Note also that for some kinematical invariants we will need labels 



from different periods. For example S31 = (1 + 3)' 



x 



35 • 



Moreover one kinematical 



invariant can have different representations from x's from different periods. For example 
the following identity holds: s 3i = (1 + 3) 2 = x^ 5 = x\ 



u 35 — ^02- 
In such notations we have for the coefficients before -B 2 " 16 

t»2...s- 



1 ,»!...«' 



grals: 



rp2m 

S2...n,q 2 



scalar inte- 



c 



Ime 
S2...S- 



\2me 
S2...S 



-l,si. 



■l,«l.. 



iytr. 
Z n 



treeMHV 



2ine 



2 S2 - 



rp rf rp nn 

• h 2s dj \s+\ dj \s- h 2s+\i 



(3.63) 



C 



2m 
S2...1 



7tree,MHV 



,2m 



-A 2 

9 S2...n,lT 
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Figure 5: Diagrammatical representation of the R rst . 



A 2m 2 



X 2n+1 X ln+l m 



(3.64) 



For example, the finite part of three point MHV form factor can be written in terms of 
dual variables as: 



/4 



fin 



1+p + p 2 ) Li 2 1 



x 



ii 



X 



13 



+ -Lo ? 



x 



l:j 



X 
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71 

~6~ 



(3.65) 



here P is permutation operator which acts on the momenta or dual variables labels. Note 
that there are no periodicity conditions on indices of x dual coordinates. We will use such 
variables for the coefficients before scalar integrals also in the NMHV case which we are 
going to discuss now. 

4 NMHV form factor 

For NMHV form factors in general we have the folloving expantion in terms of scalar 
integrals: 



y(l),NMHV \ "* /i2me n2me , si2mh r>2mh , /^i3m p3m , \ "* /^i2mrp2m _i_ /~i3mrp3m , 

n / j % i % i i % / j J J J J 



perm. 



(4.66) 
(B lm scalar integrals can also appear as the limiting case of B 2me ). NMHV form factors 
should be the next-to the lowest components in 77's expansion of Z n of Grassmann degree 
8+4, so the Grassmann degree of the C coefficients are 12. This implies the following 
configurations of MHV and MHV vertexes 



MHV x MHV x MHV x MHV 3 , 
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Figure 6: Diagrammatical representation of the R, 



(2) 

rst' 



and 



NMHV x MHV x MHV3 x MHV, 



for the quadruple cut integrands which defines the coefficients before box type scalar 
integrals. The configuration involving NMHV vertex (which can be amplitude or form 
factor) can be treated recursively. To compute n point NMHV form factor at one loop one 
will need n — 1 point NMHV tree form factor (all the tree NMHV amplitudes are known 
at least in principle in the form which can be used in our computations [33j[36]). One can 
extract n — 1 point NMHV tree form factor from n — 1 point NMHV one loop form factor 
using f)3.6ip . The recursion starts with the MHV3 form factor which is also NMHV3 one. 
In the NMHV 3 case there are no contributions from NMHV x MHV x MHV 3 x MHV 3 
configuration. So one can extract from (13.611) the form of NMHV3 at tree level and then 
use it in the computations of NMHV4 from which using f)3.6ip one can extract NMHV4 
at tree level and then use it in the computations of NMHV5 ect. One can also obtain 
the form of NMHV3 at tree level using the representation of NMHV3 as MHV vertex 
in conjugated r\ variables [36]. We will use both methods as consistency check and will 
discuss the structure of MHV3 = NMHV3 form factor in the next sub section at tree and 
one loop level. 

For the triple cut integrand situation is the similar and one have to consider the 
following type of integrands: 



MHV x MHV x MHV, 



and 



NMHV x MHV x MHV 3 . 

The last case can be considered in full analogy with the quadruple cut case. The case of 
triple MHV cut is more complicated and one have to use explicit kinematical solutions of 
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However in the case of T 2m triangles at least in the case of n = 3, 4 one can take a 
short cut and fix the coefficient using the (I3.6ip without any direct computations. 

Let's now discuss the general structure of MHV x MHV x MHV x MHV3 quadruple cuts 
which correspond in general to the coefficients C 3m before B 3m integrals (the coefficients 
before B 2mh and B lm in general can be obtained using different combinations of C 3m 
coefficients and cuts involving NMHV vertexes just as in the amplitude case [36J ) . There 
are two types of configuration of vertexes which we will call 1) and 2) (see figE] and 
E]) and the special limiting case of 1) (see fig|7j). We will discuss this last case in details 
other cases can be considered in the same fashion. The coefficient C 3m 2 before 

*r+l...t— li*»...t— li<T 

B 3m 2 scalar box integral is given by the following expression: 



t-1 



r+1 



S 




r-1 



Figure 7: Diagrammatical representation of the R, 



5(i) 



rtt- 



S r + l... t -l,St...r-l,q 2 



X 



n« 4 ",? 



<%«(E-, 4 , fe \vt + 7")^ R (E-, 4 , b \'fl 



i=l 



(huy 



Hrhfoh} + ViAhr] - r)h [rk]) ^(E-l'+i Kvf + Kvf 2 ~ Kv( 



[rklfahWir] 

Kvf + Kvf 4 



(l 3 l 2 ) ...(s- I/3) 



($8/y^--l \ „A 1 \. ^A 



KvC 



(hk) . . . (r - Ik) 



(4.67) 



This expression can be simplified by using the following strategy similar to those used 
in [SB]: one has to split the third and the fourth 5 8 functions into products of two Sq r 
(or use (I2.32p prescription) and add their arguments to the arguments of firs two Sq^s. 
After that the dependence on rj^ cancels out and we obtain the Sq r 's with the total super 
momentum as their argument. The remaining delta functions can be integrated using 
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(12.30p which after some simplifications gives the following result: 



/^i3m _ iytree,MHV 

L^o . . . . .. . „2 — 6„ 



(12) . . . (nl) ( [hr] x ; 



■ '■■' ' • '■"■ " den. \[l 3 l 4 ] 

r—l s— 1 



(4.68) 

where den. is the product of all denominators in (14.67)) . The argument of the S 4 can be 
simplified using the kinematical constraints of MHV3 vertex which implies relation 

(i\X\lx) = (i\X\r)frfr, X = hh,hh, X { is arbitrary. (4.69) 

[hh\ 



So using momentum conservation conditions associated with the MHV and MHV3 vertexes 
we obtain: 



,3m 7 tree,MHv ( 12 )---( nl ) ( MM 

" M ' -"■' : " /J ' '""' den. \[hk][kk] 

'r-l s-1 



<5 4 ^2vi( i \ ( lPr...t-i\r) +y^ J Vi(i\qPt...r-i\r) 



i=t 



(4.70) 



which finally can be written as: 



1 



(~iZm _ rrtree,MHV 6(!) 2. A 

( - y s r+ i... t _i,s t ... r _i,g 2 — ^n -"Tit 9 *H-i...t-l,s«...T— 1,^ 

£(1) __ (g ~ !)^ (E[=t 1 ^(%Pr...f-l|r) + E'=r^( j |gPt-r-lk)) 

g 4 (r|p r ...t_ig|i)(r|p t .. r 5|t - 1) (r\x t .., r -iq\r) 

^*r+l...*-l,s«...r-l,g 2 = S r...i-lV..i ~ Si...r-lS r +l. ..t-1- (4-71) 

This expression can be further simplified and written in a compact form if we introduce 
dual Grassmann coordinates: 



t-l i-1 

'rt\ — 



5>i(*l=X> A *' (4-72) 



«=r 



^-y3m ytree,MHV r>( 1 ) _ A 

U s r +i...t-l,st... r -l,q 2 ~ n ft-rtt 2 s r+i...t-i,s t ...r-i,q 2 

5(1) (tt - 1)(5 4 ((^.^in+iXri |r) + (fl ra |Sin+lZfr|r)) 

^ln+lV Prt^ln+lF/Vptr^ln+ll^ ~~ 1) ( r ptr^ln+lr) 
^s r +l...t-l)**...r-l,9 a = X rt X rt+l ~ X tr X r+lf V*. ( Oj 
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Note, that the labels r, s, t when we are using standard helicity spinor notations belong 
to the corresponding momenta, and obeys periodicity conditions i + n = i, i < n. But 
when we use dual coordinates x, 9 which lives on the infinite periodic contour there are no 
periodicity conditions on the r, s, t labels any more. Throughout this paper if not men- 
tioned otherwise we will think of r, s, t as labels belonging to the corresponding momenta, 
and so obeys periodicity conditions i + n = i, i < n. 

Similar computations in the 1) and 2) cases give the results: 

'-"' S s -l... r + lS t ...r-lSt...s-l fl - t * J rstr ) L ^Ss-l...r + lSt...r-lSt...s-l-l 



2 



d(1) 



(s- ls)(t- lt)5 4 ((8 tr \x ts x rs \r) + (O rs \x ts x tr \r)) 



x 2 s (r\x rs x ts \t - l)(r\x rs x ts \t){r\x t rXts\s - 1) (r\x tr x ts \s) ' 

^■S s -i... r +i8t...r-ist...s-l = s r...s-lS r ...t ~ St...r-lS r +l...s-l = X rs X rt+l ~ X tr X r+lsi Y^-'^l 



^3m ___ ytree.MHV f>( 2 ) _ A 



-T. . .3 



St..r—lSt...r—l n ''rst r\ s r...sSt..r-lSt...r—11 

"4 



d(2) 
Jri 'rst 



(s - ls)(t - lt)<5 4 ((6 t r\x st x sr \r) + {6 rs \x st x t r\r)) 



x 2 st (r\x sr x st \t - l)(r\x sr x st \t)(r\xt r x s t\s - l)(r\x tr x st \s) ' 

^■s r ... a st..r-ist...r-i = S s. . .r—l s t. . .r ~ Sr...t—l s s...r = X sr X tr+ i — X rt X sr+ i. [4:. 10) 

Note that these expressions are very similar to those of the coefficients R rs t in the NMHV 
amplitudes [36] • The structures of R rst , R rst and R rtt are in fact identical and the only 
difference is the rearrangements of the sums in x and \9) dual coordinates which are made 
in such a way that to avoid the dependance on q and 7 axillary variables which parame- 
terize the dependence on the (super)momentum of the operator. Such rearrangements are 
always possible due to the total (super)momentum conservation conditions in z^ ee ' MHV . 
In the case of amplitudes there are large sets of relations between different combina- 
tions of R rst [S3] • For example 

Rr,r+2,t — Rr+2,t,r+l- (4.76) 

For the case of form factors one can show that for n = 4 such relation gives 

tt r ,r+2,t = -*V+2,t,r+l- (4.77) 

To see this one has to consider relation R r>r +2,t = R r +2,t,r+i for n = 6 and then put 
q^ = Xr] a ,qQ=Xr] a and q§ = q$ = 0. The clockwise cyclic order of external legs is 
implied. Such relations in the case of amplitudes can be easily proved using momentum 
twistor representation [49|. It is interesting to note that obtained here R^ st and R r 



L rst aliu - ± L rtt 

coefficients can be rewritten in momentum twistor notations as well, and are equal to 
special cases of [abcde] momentum twistor invariants [1H] . We are going to discuss this in 
more details in separate publication. 
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In the NMHV computations one will also encounter the T 3m 2 three mass tri- 

angles. The coefficients C 3m 2 before such triangles should be fixed by the triple 

MHV cuts. In such case the explicit solutions for A^ and lf a should be used, so it is 
problematic to obtain representation for such coefficients in terms of only Aj, \ spinors, 
which corresponds to external momenta. For the integrand of such cut one has: 






(Ilk) ■ ■ ■ {sh) 



(4.78) 



(hl 3 ) ...{s + l/i) 
After the integration over Grassamann variables one obtains: 

7 MHV,tree W(^ + *)g EU^Wi^) + EL + 1 »ft W M) ,. ?q) 

(rhHWiiisHs + ihHhkHhtHhh) 2 ' { ' ' 

At this moment one has to use explicit solutions for A; r After substitution of these 




s+1 



Figure 8: Diagrammatical representation of triple MHV cut. 
solutions one can take t — ¥ oo limit and obtain: 



/-r3m 

s r ... s ,s s+ i...t,q 2 



K 



rst 



yMHV,tree ■ 



n 



■7v. r .s*A 



rst L -±s r ...s,s s +l...t,q ) 



±s 



7 (K\ 



K\K\ V 7 



— -1 



-3 



(rt){ss+ 1) 



x 



5 4 tf?/ 7 j></#) + {Klh ~ 1) £ **<*& 



z=s+l 



,2 2 



^s r ... s ,s s+ i...t,g 2 — 9 — ^ln+l- 



(4.80) 
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Figure 9: Contributing cuts for the box B lm and triangle T 2m scalar integrals coefficients 
to NMHV3 form factor. Permutations of external momenta are not shown. 

Here K 1 = p s +i...t = 3Vfi...t-i> K 2 = q = a?i n +i and K\ - are the massless projections of 
one massive leg in the direction of another masslessly projected leg (see appendix). 

Let's now discuss the beginning of the recursive procedure discussed above for the 
NMHV tree level form factors which one will need in general for computation of quadruple 
cuts. We will discuss one step of this procedure and obtain answers for n = 3 and n = 4 
point NMHV form factors at tree and one loop level. 

4.1 3 point NMHV form factor at tree and one loop level 

Let's consider the representation for MHV3 form factor at tree level in terms of 77 variables 



{ W 3 \T%\0) MHV = Ff HV ({\,\fj}) 



Fgr (EjU Kvt 

[12] [23] [31] 



where 



y(°) d6 = (7=^6)2 



kab\2 



Jab 



kab\2 



(4.81) 



(4.82) 



We can perform the Fourier transformation from 77 to 77's [3S] and write representation 
pMHV f or Mjjy 3 form factor at tree level in terms of 77's as: 

Ff^dXXv}) = j\{d% expivffjiW^axXf}})- (4-83) 



This integral can be evaluated exactly and for Z^ IHV which is connected with F^f HV as 

-ptreeMHV = exp ^a g «)_prtree,MHV 

z treeMHV = f^tree,MHV^ (4>g4) 



2(3 



we can obtain, noticing that MHV3 form factor is also NMHV3 one: 
zr ,MHV = zr , NM nv = SU J2 XX + 7 «) f[ 5^-^^^^- (4-85) 

This expression is explicitly cyclically invariant. Using identities: 

(fli2|xi4X 3 4|2> + (#24|xi 4 x 12 |2> = (12) (23) (77! [23] + 7? 2 [31] + 77 3 [12]), (4.86) 



and 



3 
4 



X 14 



W 2 (rjt) = S' GR (J2 X ^ 2 mm +Vm + %1l2]), (4.87) 



., rytree,NMHV 

one can write Z 3 as: 

rytree,NMHV rytree,MHV 5(1) /, oc ,\ 

Z 3 = Z 3 H 2U , (4.88) 

where i? 211 * s given by f)4.7ip for n = 3. Note also that in this case (n = 3) the following 
identity holds: i? 211 = -R 322 — -^133 which is consequence of the cyclical symmetry of the 
initial expression for NMHV three point form factor at tree level. Using this we can write 
Z 3 ree ' in manifestly cyclically invariant form using R rtt coefficients: 

z tree,NMHV = ^IHV^ree^ + p + ffi)^ ( 4 . 89 ) 

o 

where P is permutation operator. We can now use Z 3 ree ' in the computations of 

Z4 , but before that let's perform one loop computation also for the Z 3 ' ' , 

which is essentially trivial. The contributing scalar integrals are the same as in the MHV3 
case: B lm and T 2m . Moreover because NMHV3 = MHV3 the ratio of the one loop 
correction over tree result will be the same as in MHV3. The coefficients C lm before B lm 
integrals are fixed by the quadruple cuts of the type (14.711) . For example for B lm (l, 2, 3|g 2 ) 
integral we obtain: 

r \m _ 7 tT-ee,MifV 1 5(1) 2 2 
°q 2 — ^3 9 211 24 13 ' 

The computation of the C 2m coefficients before scalar integrals T 2m is more complicated. 
The integrand of corresponding cut is given by the triple product of MHV vertexes, and 
there are not enough kinematical conditions to express lf a momenta and A/ a spinors in 
terms of the external momenta and spinors associated with them not using exact form of 
±5* solutions. However there is only one type of such coefficients and one can fix their 
value by requiring the (I3.6ip condition must be satisfied. We will use this approach from 
now on. 
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Figure 10: List of all contributing cyclically inequivalent scalar integrals to NMHV4 



One can see that in the case under consideration the (13.611) condition satisfied if and 
only if (let's consider coefficient before T 2m (l\q 2 , (2 + 3) 2 ) integral): 



U S23,q 2 



7 tree,MHV 1 5(1) / 2 _ 2 
^3 - fl 21ll x 24 x 



UJ 



Coefficients before other integrals can be obtained by action of permutation operator 
P. Combining all contributions together we can arrange the final result as: 



y(l),NMHV iytree,NMHV 



i# lm (l, 2, 3|g 2 ) + i# lm (l, 3, 2\q 2 ) + ±B lm (2, 1, 3|g 2 ) 



+ T 2m (l|g 2 , (2 + 3) 2 ) + T 2m (2|g 2 , (1 + 3) 2 ) 
+ T 2m (3|g 2 ,(l + 2) 2 ). 



(4.90) 



Substituting the expansions in e of B lm and T 2m functions we obtain the (I3.65P result. 



4.2 4 point NMHV form factor at one loop 

Now we are ready for the computation of NMHV4 form factor at one loop. The contribut- 
ing scalar integrals in this case are 11 ! B 2mh , Q lm ^ r p 2m anc | j"3m rj^g g rs ^ ^ ree types of 
integrals are IR divergent, while the last one (tree mass triangle) is IR finite. Let us 
remind the reader that there are no other contributing scalar triangle integrals for the 
reasons discussed in the previous section. Let's concentrate first on the IR divergent part 



11 



The B 2me scalar boxes are absent by the same reasons as in the 6 point NMHV amplitude 
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of the answer. We will label the coefficients before cyclically inequivalent type of integrals 
as: 

/~t2mh /~t2mh film ri2m (a q-t \ 

°« 2 (l+2) 2 ' '-V, (4+3) 2 > ( - / (<7-4) 2 ' ( -V,(l+2+3) 2 ' l^.yij 



y g 2 ,(l+2) 2 ) ^ ? 2 ,(4+3) 2 ' ^(q-if'^q 2 , (1+2+3) 

Tim l^\ nirn ri 

V,(l+2) 2 

% 2 ,(l+2) 2 



For Cg™({ i 2)2 we have two types of cuts contributing (see figJTTT): 



C fl 2 2 t 1+2)2 =A) + 5). (4.92) 



r^nA - ytree,MHV L i 6(1) , p(2) \ 2 2 /. q o\ 

°g 2 ,(l+2) 2 — ^4 o ^311 + n 4XZ) x 35 X lA- ^4.y<3J 



These cuts give us: 

■ 1 , a 

V,(l+2) 2 _ ^4 2 

For Ci™ 4 y 2 we have two types of contributions: 

C\™, ? =C) + D). (4.94) 

These cuts give us: 

film _ 7 tree,MHV 1 /6(l) , d(1) \_.2 2 /, q rN 

°(g-4) 2 — ^4 q^ 144+ 241 ' 24 13 ' ^.yOJ 

Note that C) cut contains vertex which is NMHV 3-point form factor. To evaluate such 
cuts we use the same technique as in |36j: we substitute explicit expression for NMHV 
3-point form factor and then use the kinematical constraints attached to the neighboring 
MHV3 vertexes which imply A; 2 ~ Ai and A; 3 ~ A3. Using the cyclical symmetry of 
NMHV 3-point form factor we can obtain the following relations for R^: 

6(1) _ 6(1) 6(1) _ 6(1) ( A QK\ 

-"-144 ~~ -"-311' -"-244 ~~ -"-211 • \*.VV) 

For (7^(4+3)2 similar to C^f 1+2)2 we have: 

r 2mh _ 7 tree,MHV ± / 6(1) , p(2) \J2 2 ( A Q7\ 

L 'g 2 ,(4+3) 2 - Z 4 9^ 244 142 ' ) 13 25 ' l 4 - y 

For C 2 Tii , 2+3) 2 we have two contributing cuts. Cut involving NMHV 5 point tree 
amplitude can be evaluated using explicit expression for NMHV 5 point tree amplitude 
and kinematical constraints associated with the neighboring MHV3 vertex, which imply 
\i 2 ~ Ai. We will adjust the value of the other triple MHV cut using the A3 . 6 1 [) just as in 
the previous NMHV3 case. Finally we have: 

r T2m _ y tree,MHV*-/nQ.) , d( 2 ) \( „2 2 \ (a qq\ 

O g 2 i(1+ 2 +3 )2 - Z 4 -( K 311 + K 413 )[x u - x 15 ). (4.98) 

One can notice that CrfFr^+g) 2 — ^^,? 2 ™(i+2+3) 2 ' so we can WT ^ e an contributions from 
the IR divergent triangle integrals in the following form: 

C 9 2m (l+2+3) 2 ^ 2 7l+2+3) 2 + ^ 2 ™4+2+3) 2 ^ 2 74+2+3) 2 + P erm - = 2C g2 7;i + 2+3) 2 '^ 2 "(l+2+3) 2 + perm. 

(4.99) 
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So, finally gathering all IR divergent contributions together we obtain 



r NMHV,(l) , 7 tree,MHV 

'4 / Z 4 



1 ,~ 



(1) 



? (2) \ R 2m/i 



IB 



(R£ 1 + R% 3 )B^(3, 4|^,(1 + 2) 



IR 



IR 



+ ^( J RS 1 4 ) 4 + 4 1 4 ) i)^ lm (l,2,3|(g-4) 2 ; 
+ ^?44 + fl&^a, 2|g 2 , (3 + 4) 2 ) 
+ |(flft } i + 4i ) 3 )2T 2m (4|g 2 , (1 + 2 + 3) 2 ) 



IR 



IR 



+ perm. 



Using equation 

which is consequence of R r , r +2,t = R r +2,t,r+i 



,(2) 
L 413 



iX/MQ -ft. 



(1) 

211 



(4.100) 
(4.101) 



relation for amplitudes, eg. (I4.96|) and 



5(i) _ 6(i) _ ip 6(i) 

-n, 2 44 — -fl-211 — Jr - rL 144? 



also noticing that R^ 2 = P-R241 we can write: 



Y NMHV,(1) ,rvtree,MHV 
^4 7^4 



1,*< 



;(4 1 i ) i + 4 1 4 ) i){^ 2m/l (3,4|g 2 ,(l + 2) 2 ) 
+ B 2mh (A, l\q 2 , (3 + 2) 2 ) + B lm (l, 2, 3\{q - 4) 2 ) 



+ 2T 2m (4|g 2 ,(l + 2 + 3) 2 )} 



IR 



+ perm. 



(4.102) 



Using explicit expressions for the IR divergent parts of integrals 



£ 2mh (3,4|g 2 ,(l + 2) 2 ) 

£ 2mft (4,l|g 2 ,(3 + 2) 2 ) 

£ lm (l,2,3|(g-4) 2 ) 

T 2m (4|g 2 ,(l + 2 + 3) 2 ) 



/' 



IR 



IB 



IB 



IB 



e 2 


^ 


e 2 


^ 


e 2 


^ 



(/& + 2(4 - g)* - g* - *y 

(^ + 2(4-^-^-^) 
(2 S | 2 + 2 S21 -2(4-g) 2e ), 
(q 2e - (4 - q) 2 ') , 



(4.103) 



one can see that IR divergent part of the NMHV four point form factor indeed has the 
form: 



7 NMHV,{1) irrtree,MHV 



/z\ 



IB 



5(1 + P + P 2 +0(^1 +i^)^ 1 fS " +l 



j=l 



e 2 \ fi 2 



(4.104) 
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From this we conclude that 



z NMHV,tree = ^HV.tree^ + p + p 2 + p3) (£« + ^ (4 _ 1Q5) 



Let's now consider the contribution of IR finite T 3m ((l+2) 2 , (3+4) 2 , q 2 ) scalar integral. 

'(1+2)2,(3+4)2 ,<j2 



For the corresponding coefficient C?^ 2 s 2 (3.4)2 „2 using previous results we immediately 



obtain: 

^(l+2)2,(3+4)2, g 2 

The finite part of four point NMHV form factor then can be written in the following form: 



q 3 i+2) 2 ,f3+4) 2 ., 2 = Z MHV ^K luq \ (4-106) 



z NMHV,(l) /z tree M HV = ^ + p + p 2 + p3) I ( £(1) + jjWj ^ 

fin 2 



+ (i + p + r + F% 4 iy 4 . (4.107) 



Where (Log(:r) = L(x)) 

2 



V 4 = -2^fLi 2 fl-^+Li 2 fl-^)-4Li 2 fl-^ 



9 \ / 9 \ / 9 \ \ / 9 , 

-y \ I rtr*^ \ / rf^ \ \ I ry^ \ IT 



i=1 V V x 14 / V x j+2J+4/ V x i+2i+4/ / V x 24/ ° 

(4.108) 

while Wi is given by Davydychev function [50], which in our case has the form: 
W 4 = ^(2U 2 (-xR)+2U 2 (-yR) + L(xR)L(yR) + L(^L(^±^+^ 



1 / 2 D n(i „ „. i dW ™ ^13 „, X 35 



Q = {{l-x-y) 2 -Axyy'\ R = 2{l-x-y + R)-\ x = ^f,y=^f. {Aim) 

x 15 x 15 

This is the end of computation of the four point NMHV form factor. Using results 
obtained in previous section for MHV sector one can in principle define analog of R£ 
ratio function as in the amplitude case |36J . One also can use (I4.105P as an input in the 
computation of the five point NMHV form factor. The structure of such computation 
should be essentially similar to the four point case. The IR divergent part of the answer 
will be determined by the IR divergent parts of box and two mass triangle scalar integrals, 
while in the finite part there will be additional contributions from IR finite three mass 
integrals. 

It is interesting to note that the number of R^ st coefficients in tree level 3 and 4 
point form factors exactly matches the number of diagrams which one will have for the 
corresponding super Wilson loop [25] if one will use the prescriptions of [13] about selection 
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of diagrams. It will be interesting to clarify by explicit computations the status of the 
Wilson loop/form factor duality in the none-MHV case. 

In [171 EH] it was observed that the following relation between form factors and am- 
plitudes likely holds 

z n ({\, a, n}, o, {o}) = t[T? HV ){{\ a, 77}, o, o) = g ^M\\v}_ (4110) 

This relation was verified at tree and one loop level for the MHV sector. While for the 
component answer for the 4 point NMHV form factor the limit q — y was in general 
singular. If one considers the coefficients before corresponding scalar integrals in the 
case of five point form factor at one loop one can observe the following: in the limit q —y 
0, 7^ — y the coefficients before triangle integrals vanish and so do most of the coefficients 
before box scalar integrals. The only coefficients that survive are the coefficients before 
B lm box integrals which in this limit reproduce the coefficients for the five point NMHV 
amplitude. It is likely the general pattern for the n point case. One may think that 
relation (I4.110p holds for such numbers of external legs for which the objects on both 
sides of the equality exists. I.e. equality should not hold for example for the NMHV 
sector for n = 4 because there are non = 4 NMHV amplitude. 

5 Conclusions 

In this paper the systematic study of form factors in the M = 4 SYM theory in on-shell 
momentum superspace formalism is performed for the MHV and NMHV sectors at one 
loop order of PT by means of generalized unitarity technique. The use of M = 4 covariant 
methods allows us to obtain answers for any type of operator from stress tensor multiplet 
and arbitrary external states in these sectors. The explicit answers for the 3 and 4 point 
NMHV form factors were obtained as well as the n point situation was discussed. As the 
byproduct of this investigation the representation for 3 and 4 point form factors that does 
not depend on any kind of "reference spinors" at tree level were obtained. 

The application of the generalized unitarity methods to form factors clarifies several 
issues: the structure of the basis of scalar integrals at one loop and the relations between 
form factors with operator insertion with zero momenta and the amplitudes. 

The recent studies of the structure of the amplitudes and their relations to Wilson 
loops in M = 4 SYM led to the formulation of the equation which at least in principle 
should define full S-matrix of the theory [2H1 HD] at any value of the coupling constant. The 
conjectured derivation of this equation is based on the amplitudes/Wilson loops duality. 
There is also the conjecture [13] that similar duality between form factors and Wilson 
loops also holds. It is interesting to investigate whether such duality survives for the 
NMHV and other sectors in some form and if it is possible to formulate similar equation 
for the form factors. The results obtained in this paper may be considered as starting 
point in such investigation. 
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A J\f = 4 harmonic superspace 

We discuss here the reformulation of f )2.19p in the M = 4 harmonic superspace. Our dis- 
cussion is based mostly on section 3 of [32J . The M = 4 harmonic superspace is obtained 
by adding additional bosonic coordinates (harmonic variables) to the M = 4 coordinate 
superspace or on-shell momentum superspace. These additional bosonic coordinates pa- 
rameterize the coset 

SU(A) 

SU(2) x SU{2)> x U(l) 
and carry the £77(4) index A, two copies of SU(2) indices a, a and £7(1) charge ± 

«, v). 

Using these variables one presents all the Grassmannian objects with SU(4)r indices. 
For example, for Grassmannian coordinates in the original M = 4 coordinate superspace 

0t a = <0t K a = u~ A a 6t (A.Hl) 

and in the opposite direction 

0t = + a a u A +a + 6~ a u A +a . (A.112) 

The same can be done with supercharges etc.. Note that harmonic variable projection 
leaves helicity properties of the objects unmodified. Also, similar projections can be per- 
formed for Grassmannian coordinates r\ A and supercharges q A , q^A of on-shell momentum 
superspace. 

So the M = 4 harmonic superspace is parameterized with the following set of coordi- 
nates 

N = A harmonic superspace = {x a& , 6+ a , 6~ & , 9+ a , 9^, u] 

or {\ a ,\a, V +a ,V' d , u}. (A.113) 

Using u harmonic variables one can project the W AB superfield as 

W AB _^ W AB u +a u +b = e ab w ++^ 
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where e ab is an £77(2) totally antisymmetric tensor and the Grassmannian analyticity 
conditions 1321 such thato 



Dt aW ++ = 0, D% a W ++ = 0. 

Thus, the superfield W ++ contains the dependence on half of the Grassmannian vari- 
ables #'s and #'s. 

W ++ = W ++ (x } e+\6T\ u ), 

Performing the expansion of W ++ in u all the projections like (I2.8P in SU(4)r covari- 
ant fashion can be obtained. This is the main purpose of introduction of the harmonic 
superspace. The component expansion of W ++ in #'s and #'s can be found in [32]. The 
lowest component of the W ++ expansion is 

W ++ (x,0,0,u) = 0++, 0++ = ^e ab u+ a u+ b <t> AB , 

where according to [32J 

Q1 a <f> ++ = 0. (A.114) 

Using this condition and the translation invariance we can write the expression for 
MHV "super state -super form factor" at tree level in harmonic superspace |14j : 

Z^ MHV {{\ A, V }, q, u, 7 +") = ^ +4( ( f 2 + 7 " ) ( ^ (C) , (A.115) 

where S ±4 is the Grassmannian delta function; <5 ±4 are defined as 

n 2 

5 ±4 (g^)=E II Mvt a/ %? /h - (A.116) 

i J =1 a/a,b/b=l 

We can also define 5 ±2 as usual Grassmann delta functions: 

2 

5 ±2 {X a/h ) = ] [ X a/d . (A.117) 

a/a, b/b=l 

The obtained expression for form factor looks just like (I2.19p . but now both the Grass- 
mannian delta functions 5 ±A are SU(4)r covariant. One can write also the MHV part 
of a superamplitude in a similar manner. Projecting the condition of superamplitude 
invariance under q£ supersymmetry transformations we have 

q AJtree,MHV = Q ^ {g +a + q - a) £ree,MHV = ^ 



12 Strictly speaking this is true only in the free theory (g = 0), in the interacting theory one has to 
replace D^,D^ by their gauge covariant analogs, which contain superconnection, but the final result is 
the same [32]. 
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and taking into account that the helicity properties of projected supercharges are not 
modified we get 

Atree,MHV " \Qa)° {QaJ f\ 110 \ 

An = (12)... (nl) • (A ' 118) 

Now both ^^■eei-W'-ff^ anc [ ztree,MHV are gjj(^ R invariant and one can use them in uni- 
tarity based computations, where Grassmann integration (super summation) should be 
performed separately for d +2 r] and d~ 2 r\. 

We see now that all results obtained in previous sections can be simply generalized to 
SU(4)r covariant harmonic superspace version. We have to replace common z^ e 
prefactor by its harmonic superspace generalization and replace all 5 A {X^) functions in 
Rf st by combination of 8' 2 (X- & )S +2 (X+ a ). 

B Kinematical solutions 

This chapter is based on [37]. For the case of triple cuts for the cut momenta lf a and 
associated with them spinors one can find explicit expressions in terms of external mo- 
menta data and t parameter which is the remaining of the loop integration. First one can 
define massless projections of combinations of external momenta: 

, aa Kf ~ KfKfh , a6l Kf - KfK 2 h 
Kl ~ 1-KlKlh 2 ' K2 ~ \-K\Klh 2 • [ 9j 

Here 7 = (Ki, K 2 )±((Ki, K 2 ) 2 — KfK 2 ) 1 ^ 2 which corresponds to two possible kinematical 
solutions ±S. In general K 2 7^ 0, j = 1...3. 

Using these massless projections one can define corresponding spinors. The notations 
of [37] X a , = {K\\ = {K\-\ and A° b = [K\\ = {K\ + \ are used in the sense that: 

i i 

(ij) = ( r \j + ) =u-(ki)u + (kj) } [ij] = (i + \j~) =u + (k i )u^(k j ), 
where u± are four component Weyl spinors. Then we have: 

(l-\ = t(K\-\+a a (K b 2 -\, (l+\ = ^(K\ + \ + a a (K b +\, (B.120) 

while for lf a = a" a /f one can write: 

If = a i2 KY + aaKf + ^M 7 ^-) + < ^{k\-\ 1 »\K\-). (B.121) 

Here the explicit expressions in terms of K 2 and 7 for ctij can be found in appendix A of 



Kl^-K 2 ) _ K^-Kl 



«01 — — o TPTTPTi a 02 



! 2 -K 2 Kr "- ^-iqiq 
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"ll = "01 , "l2 = "02 - 1, 

7 

K 2 

"21 = "01 - 1, "22 = "02 -■ (B.122) 

7 
Using these expressions one has the following relations for the spinor products of (klj): 

(hh) = -t(l-^{K\K\), 

tK 2 
(hh) = ~^{K\K\), 

7 

(l 2 l 3 ) = t{K\K\). (B.123) 

In the case when K 2 = fore some j (T 2m scalar integral) the explicit solution is 
different. If the massless leg is attached to the MHV vertex then the solution takes the 
form (here we assume that K\ = 0): 



t /K+ \ K K 2 [ X K\] 

WM ; 2 lKl + - {x l{k ^WM { 2l 



[K 2 K\] X 2l Vl1 [K.KW 

1 



<TI = TTT^i^lKi, (lt\ = -t^^(Kt\ ' (K\\ 

<'2"l = T7^( K l\ K 3 + 77^( K 2\Ki, (ll\ = (ltl (B.124) 



where (x\ is arbitrary spinor. Using these solutions we have the following relevant for our 
computations results: 

(hh) = -^-(xKl), 

7 

(hh) = - — (XK\), 

7 

(hh) = [ -^3(hK 2 ). (B.125) 

\X K il 

C Scalar integrals 

For the Box type integral: 



d D l 1 

'KiXiXIXl - J ( 27r )^ p(K[ + /)2(^ +K2 + /)2 (/ _ Ka) 2 > 
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(C.126) 



where D = 4 — 2e we define dimensionless function B K 2 K i K 2 K 2 as: 

&KIKIKIKI = {iK D/2 rrY l A(2tt) d B k1k1k1k! , (C.127) 

where 

r(l-e)r(l + e) 
r r = 



r(i - 2c) ' 

while for A for B 3m , B 2mh } B 2me and B lm box scalar integrals we have (s^ = (iTj + Kj) 2 ) 



a 3m 


- S12S23 " 


- K 2 K 2 


A2m/i 


= S12S23, 




A2me 


= S12S23 ' 


- K 2 K 2 


A lm = 


= Si 2 S 2 3- 





(C.128) 

For the B functions we have [31] (note that we rearranged IR divergent part of B 2mh in 
comparison with 



e \ s 12/ V S 12/ 



2L,, ( 1 _ S] _ 2Li 2 (l-£±\- Log 2 ( S -^) CHr) 



s 23/ \ s 23/ \ S 23 



B 2 ™(l,3\KiKJ) = ^(2si 2 + 2 S tl-2K?-2K?)-2Li 2 (l-^\ 

-2Li 2 (l-%)- 2Li 2 (l-Z±\- 2Li 2 (l-*l 



■S23 / V s 12 / V s 23 

■Sl2 s 23 / \ S 23/ 



B lm (l,2,3\K%) = ^(2sl 2 + 2s^-2K?)-2U 2 (l-^-)-2uJl-^- 

£ 2 \ S 12/ V S 23 



-Log 2 ( S f) -y + 0(e). (C.129) 



For triangle scalar integrals 

f d D l 1 

T Kt,KlKl = J {27[) D HKl + m _ K ^ r ( C ' 13 °) 

we have similar definitions: 

TkIkIkI = {i<K D/2 r T y l A(27r) D T KlKlKl (C.131) 
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While for A coefficients for the relevant for our discussion cases (T 2m and T 3m scalar 
triangles with q 2 massive leg) we have: 



3m . _ 2 

. 2m 



A 2m = K *_ q \ (C.132) 



r 2m (l\K!,q 2 ) = ^(Kt-f). (C.133) 

The T 3m triangle is IR finite and the answer for it is given in terms of Davydychev function 

T im (KlKlq 2 ) = r 3m (K 2 / q 2 ,K 2 / q 2 y. 

r3m = 2U 2 (-xR) + 2Li 2 (-yR) + L og(xR)Log(yR) + Log (g) Log (jgg) 




3Q' 



7T 2 



g = {{l-x-yf-Axy) l/ \ R = 2(l-x-y + R)-\ x=% y=%. 



(C.134) 
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Figure 11: List of contributing quadruple and all possible triple cuts for NMHV four point 
form factor. B), G), D), F) triple cuts give vanishing results for corresponding coefficients. 
Permutations of external momenta are not shown. 
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